Predicted modulated differential rates for direct WIMP searches at low 

energy transfers 

J. D. Vergados * 
CERN, Theory Division, 
CH 1211, Geneva 23, Switzerland 
and 

Theoretical Physics Division, 

University of loannina, 
loannina, Gr 451 10, Greece 
(Dated: September 2, 2011) 

The differential event rate for direct detection of dark matter, both the time averaged and the 
modulated one due to the motion of the Earth, are discussed. The calculations focus on relatively 
light cold dark matter candidates (WIMP) and low energy transfers. It is shown that for sufficiently 
light WIMPs the extraction of relatively large nucleon cross sections is possible. Furthermore for 
some WIMP masses the modulation amplitude may change sign, meaning that, in such a case, the 
maximum rate may occur six months later than naively expected. This effect can be exploited to 
yield information about the mass of the dark matter candidate, if and when the observation of the 
modulation of the event rate is established. 
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INTRODUCTION 

The combined MAXIMA-1 Q, BOOMERANG 0, DASI and COBE/DMR Cosmic Microwave 
Background (CMB) observations Q imply that the Universe is flat @ and that most of the matter 
in the Universe is Dark [6|, i.e. exotic. These results have been confirmed and improved by the 
recent WMAP data 0| • Combining the the data of these quite precise experiments one finds: 

= 0.0456 ±0.0015, f^cDM = 0.228 ±0.013, 17a = 0.726 ± 0.015. 

Since any "invisible" non exotic component cannot possibly exceed 40% of the above Qcum jll, 
exotic (non baryonic) matter is required and there is room for cold dark matter candidates or 
WIMPs (Weakly Interacting Massive Particles). 

Even though there exists firm indirect evidence for a halo of dark matter in galaxies from the 
observed rotational curves, see e.g the review it is essential to directly detect such matter. Until 
dark matter is actually detected, we shall not be able to exclude the possibility that the rotation 
curves result from a modification of the laws of nature as we currently view them. This makes it 
imperative that we invest a maximum effort in attempting to directly detect dark matter in the 
laboratory. Furthermore such a direct detection will also unravel the nature of the constituents of 
dark matter. The possibility of such detection, however, depends on the nature of the dark matter 
constituents and their interactions. 

Since the WIMP's are expected to be extremely non relativistic, with average kinetic energy 
(T) « 50 keV(mwiMp/100 GeV), they are not likely to excite the nucleus, even if they are quite 
massive ttiwimp > 100 GeV. So they can be directly detected mainly via the recoiling of a nucleus 
(A,Z) in elastic scattering. The event rate for such a process can be computed from the following 
ingredients: i) An effective Lagrangian at the elementary particle (quark) level obtained in the 



* vergados@uoi.gr 



2 



framework of the prevailing particle theory. In supersymmetry the dark matter candidate is the 



LSP (Lightest Supersymmetric Particle) [10l - ll6l | . In this case the effective Lagrangian is constructed 



as described, e.g., in Refs. [lol - fist . ii) A well defined procedure for transforming the amplitude thus 
obtained, using the previous effective Lagrangian, from the quark to the nucleon level. To achieve 
this one needs a quark model for the nucleon, see e.g. This step is particularly important in 



supersymmetry or other models dominated by a scalar interaction (intermediate Higgs etc), since, 
then, the elementary amplitude becomes proportional to the quark mass and the content of the 
nucleon in quarks other than u and d becomes very important, iii) knowledge of the relevant nuclear 



matrix elements [22, l23|, obtained with as reliable as possible many body nuclear wave functions, 
iv) knowledge of the WIMP density in our vicinity and its velocity distribution. 

From steps i) and ii) one obtains the nucleon cross sections. These can also be extracted from the 
data of event rates, if and when such data become available. From limits on the event rates, one can 
obtain exclusion plots on the nucleon cross sections as functions of the WIMP mass. The extracted 
cross sections depend, of course, on inputs from steps iii)-iv). 



In the standard nuclear recoil experiments, first proposed more than 30 years ago [2J], one has to 
face the problem that the reaction of interest does not have a characteristic feature to distinguish 
it from the background. So for the expected low counting rates the background is a formidable 
problem. Some special features of the WIMP-nuclear interaction can be exploited to reduce the 
background problems. Such are: i) the modulation effect: this yields a periodic signal due to the 



motion of the earth around the sun. Unfortunately this effect, also proposed a long time ago [25 1 



mateiy 

and subsequently studied by many authors 2^ 11], is small and becomes even smaller than 2% 
due to cancelations arising from nuclear physics effects, ii) backward-forward asymmetry expected 
in directional experiments, i.e. experiments in which the direction of the recoiling nucleus is also 
observed. Such an asymmetry has also been predicted a long time ago [s^], but it has not been 
exploited, since such experiments have been considered very difficult to perform, but they now 
appear to be feasible 34-4^. iii) transitions to excited states: in this case one need not measure 



nuclear recoils, but the de-excitation 7 rays. This can happen only in very special cases since the 
average WIMP energy is too low to excite the nucleus. It has, however, been found that in the 
special case of the target ^^''l such a process is feasible [i^ with branching ratios around 5%, (iv) 
detection of electrons produced during the WIMP-nucleus collision: it turns out, however, that this 
production peaks at very low energies. So only gaseous TPC detectors can reach the desired level 
of 100 eV. In such a case the number of electrons detected may exceed the number of recoils for 
a target with high Z 47], 48| , v) detection of hard X-rays produced when the inner shell holes are 
filled: it has been found 49] that in the previous mechanism inner shell electrons can be ejected. 
These holes can be filled by the Auger process or X-ray emission. 

In connection with nuclear structure aspects, in a series of calculations, e.g. in [sol - lsil l and 
references there in, it has been shown that for the coherent contribution, due to the scalar interaction, 
the inclusion of the nuclear form factor is important, especially in the case of relatively heavy targets. 
They also showed that the nuclear spin cross sections are characterized by a single, i.e. essentially 
isospin independent, structure function and two static spin values, one for the proton and one for 
the neutron, which depend on the target. 

As we have already mentioned an essential ingredient in direct WIMP detection is the WIMP 
density in our vicinity and, especially, the WIMP velocity distribution. Some of the calculations have 
considered various forms of phenomenological non symmetric velocity distributions 32, 33., ■38l . Isoj 
and some of them even more exotic dark matter fiows like the late infall of dark matter into the 
galaxy, i.e caustic rings 53-53], dark matter orbiting the Sun [ZTl and Sagittarius dark matter |58(. 

In addition to computing the time averaged rates, these calculations studied the modulation effect. 
They showed that in the standard recoil experiments the modulation amplitude in the total rate may 
change sign for large reduced mass, i.e. heavy WIMPS and large A. In directional experiments, in 
addition to the expected asymmetry mentioned above, the modulation exhibits two very interesting 
patterns i) its magnitude in certain directions can be very large and ii) the location of the maximum 
and minimum depends on the direction of observation. 

In the present paper we will expand the above calculations and study the differential event rates, 
both time averaged and modulated, in the region of low energy transfers, as in the DAMA experiment 
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[5 3 . |60[, focusing our attention on relatively light WIMPS |6ll - l63l |. Such light WIMPs can be 
accommodated in some SUSY models [g^l- We will focus here on the standard Maxwell-Boltzmann 
M-B) distribution for the WIMPs of our galaxy and we will not be concerned with other distributions 
si-lil, even though some of them may affect the modulation. The latter will be studied elsewhere. 
In such a context we will explicitly show that the modulation amplitude, entering both the differential 
and the total rates, changes sign for certain WIMP masses. As a result such an effect, if and when 
the needed data become available, may be exploited to infer the WIMP mass. 



THE FORMALISM FOR THE WIMP-NUCLEUS DIFFERENTIAL EVENT RATE 

This formalism adopted in this work is well known (see e.g. the recent reviews [isl l69j). So we 
will briefly discuss its essential elements here. The differential event rate can be cast in the form: 



dR I dRo , dH 



A = —TFriA 



dQ' dQ dQ 



Acosa (1) 



dQ Aijip XfJ-pJ Qo{A) du 

dH Px mt f PrV I n .2 1 dh 



dQ AvTLp \P-pJ QoiA) du 



(2) 



with with /ir ifJ-p) the WIMP-nucleus (nucleon) reduced mass, A is the nuclear mass number and 
(T„ is the elementary WIMP-nucleon cross section, ni-^ is the WIMP mass and nit the mass of the 
target. The first term gives the time averaged rate, while the second gives the modulated amplitude. 
a is the phase of the earth (a = on June 2nd) . 
Furthermore 

^ = ^J^a^F\u)^oiaV^), £ - ^J^a^ F\u)^ liay^) (3) 

with a — {\/2firbvo)~^ , vq the velocity of the sun around the center of the galaxy and b the nuclear 
harmonic oscillator size parameter characterizing the nuclear wave function, u is the energy transfer 
Q in dimensionless units given by 

"=7TT7n ' Qo(^) = KA62]-i=40A-4/3MeV (4) 

and F{u) is the nuclear form factor. Note that the parameter a depends both on the WIMP mass, 
the target and the velocity distribution. Note also that for a given energy transfer Q the quantity 
u depends on A. 

The functions 'i/o{ay/u) and 'i>i{ay/u) can be obtained as follows: 

• One starts with a Maxwell-Boltzmann distribution in the galactic frame with a characteristic 
velocity vq equal to the suns velocity around the center of the galaxy^ 



^ Strictly speaking, since an upper cutoff is introduced to the velocity distribution, equal to the escape velocity, the 
velocity distribution should be renormalized. However the normalization integral is close to one, namely 

0rerf(yeac) - 2e~^o=<:2 j/csc uj/csc 

— Veac = , (5) 



^^"^^^^ , — ) if cay. 

VTT VO 

i.e. norm 0.9989 for yeac = 2.84 
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• one transforms to the local coordinate system: 

V 

y ^ y + + (5 (sinax — cosQ;cos7j/ + C0SQ;sin7{)s) , y — — (6) 

with Vs a unit vector in the Sun's direction of motion and 5 is the ratio of the Earth's velocity 
around the sun divided by uo- The above formula assumes that the motion of both the sun 
around the galaxy and of the Earth around the sun are uniformly circular. The exact orbits 
are, of course, more complicated jH, [t^I , but such deviations are not expected to significantly 
modify our results. 

• One integrates over the velocity integration over the angles and the result is multiplied the 
velocity y = u/uq due to the WIMP flux. 

• The result is integrated from a minimum value, which depends on the energy transfer y — a^/u, 
to a maximum y = j/osc, Vcsc = "csc/wo, Vcsc w 2.84 



The result is 



J{x) ^ 



6 cos a — 2 



erf ( —X H — 5 cos a + 1 ) + erf ( x H — (5 cos a + 1 



2 / V 2 



+erfc ( -ycsc + ^Scosa + 1 j + erfc fycsc + ^(5 cos a + 1 ) - 2 



X = a\fu (7) 



where erf(a;) and erfc(x) are the error function and its complement respectively. Furthermore since 
b = 0.135 we can expand in powers of (5 and obtain: 



with 



J{ay/u) w ^'o(aVu) + ^'l(a^/lt)cosa + ^'2(a^At)cos2a (8) 



"^aix) = ^(crf(l ~x) + eri{x + 1) + erfc(l - y^sc) + erfc(yosc + 1) - 2) (9) 
= i(5 (-erf(l - - erf(a; + 1) - erfc(l - yesc) - erfc(?;csc + 1) 
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2g-(a;-i)" 2e"(^+i)' 2e"(2'<==^<^"i)' 2e~(^'==<^+i'' 



2 (10) 



The function vE'2(a;) is small, of order and it can be ignored. If, however, the experiments, which 
attempt to measure the modulation, want to go beyond the cos a term, they should consider terms 
cos 2a rather than sin a as some of them have done. 



The functions 5*0 (2;) and 'i'i{x) characterize both the coherent and the spin induced mode [7l|. We 
should note that the function vE'i(x) changes sign at some value of x, which has implications on the 
total modulated rate, a point often missed (see Fig. [T|). Sometimes we will write the differential rate 
as: 

2 

dR . Px nit / A'r \ / 2 — /1 2 1 ( '^^ 



dQ Am 



<Tn(^) V^^A'-^(^{l + H{aV^)cosa) (11) 
\p,pj QoiA) \du ) 



In this formulation H{a^/u) gives the relative modulation amplitude (with respect to the time 
averaged one). The functions 'i^oiay/u) and ^'1(0-^^), which exhibit the general characteristics of 
the differential rates, are exhibited in Figs [5] and [31 while the function H{ay/u) is shown in Fig. 
These functions are independent of the nuclear physics. They only depend on the reduced mass and 
the velocity distribution. They are thus the same for both the coherent and the spin mode. Note 
that '^i{ay^) and, consequently, H{a^Ju) can take both positive and negative values, which affects 
the location of the maximum. 



5 





(a) (b) 
Q ^keV 



Figure 2: The function ^o(a\/u) entering the differential rate as a function of the recoil energy for 
a heavy target, e.g. ^^^I, without the form factor (a) and including the form factor (b). The solid, 
dotted, dot-dashed, dashed, long dashed and thick solid lines correspond to 5, 7, 10, 20, 50 and 100 

GeV WIMP masses. 




Figure 3: The function ^i(ay^) entering the modulated differential rate as a function of the 
recoil energy for a heavy target, e.g. ^^''l, without the form factor (a) and including the form factor 
(b). The solid, dotted, dot-dashed, dashed, long dashed and thick solid lines correspond to 5, 7, 10, 

20, 50 and 100 GeV WIMP masses. 
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Figure 4: The same as in Fig. [3] for function Hiay/u) entering the modulated differential rate as a 
function of the recoil energy for a heavy target, e.g. ^^^I. Note that this is independent of the form 
factor. The solid, dotted, dot-dashed, dashed, long dashed and thick solid lines correspond to 5, 7, 

10, 20, 50 and 100 GeV WIMP masses. 




2 3 

Q ->keV 

Figure 5: The quenching factor used in this work to transform kcV— >kcVee. 



SOME RESULTS ON DIFFERENTIAL RATES 



We will apply the above formalism in the case of Nal, a target used in the DAMA experiment 



60|. The results for the Xe target are similar [61|. The differential rates ^'^'^ ^1^' each 
component {A = 127 and A = 23) are exhibited in Fig. [7][51 Following the practice of the DAMA 
experiment we express the energy transfer is in keVee using the phenomenological quenching factor 



72| . 73| shown in Fig. [S] The nuclear form factor has been included (for the ^^^I its effect is sizable 



even for an energy transfer of 10 keV, see Fig. [6]) . The differential rate for the spin mode for low 
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Q — >keVee 

Figure 7: The differential rate as a function of tlie recoil energy for a heavy target, e.g. ^^''l 

(a) and the amplitude for the modulated differential rate ^ (b), assuming a nucleon cross section 
of lO^^pb. The solid, dotted, dot-dashed, dashed, long dashed and thick solid lines correspond to 
5, 7, 10, 20, 50 and 100 GeV WIMP masses. Note that ^ is given in absolute units. 



energy transfers is similar to those exhibited in Figs[7][51 since the spin form factors are similar. They 
are, of course, simply scaled down by A^, if one takes the spin cross section, a combination of the 
nuclear spin ME and the nucleon spin amplitudes, to be the same with the coherent nucleon cross 
section, i.e. a^j' '" = 10~^pb. For the actual spin nucleon cross sections extracted from experiment 
III and [ziQ. 



see 



The functions H{ay/u) cos a for each target component are shown in Figs l9l [T3l as a function 
of a for various low energy transfers. The corresponding quantities for the spin mode are almost 
identical. We see that for certain values of the WIMP mass the modulation amplitude changes sign. 
This may perhaps by exploited to exreact information on the WIMP mass from the data. A similar 
behavior has been found by considering various halo models and different minimum WIMP velocities 



3al33|. 

Sometimes, as is the case for the DAMA experiment, the target has many components. In such 
cases the above formalism can be applied as follows: 

— — |a -> ^^^j"77TUii M — > Ui, Xi = the fraction of the component Ai in the target (12) 



Thus we get the results shown in Figs [T2] and 11611171 The corresponding ones for the spin mode are 
not expected to be the same. 
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(a) (b) 
a — > 



Figure 9: The modulation H{a^/u) cos a with an energy transfer of 1 keVee (a) and 2 keVee (b) 
for a heavy target (I or Xe). The solid, dotted, dot-dashed, dashed, long dashed and thick solid 
lines correspond to 5, 7, 10, 20, 50 and 100 GeV WIMP masses. Note that for some wimp masses 
on June 2nd the amplitude becomes negative (location of minimum rate). Note that the 
modulation is given relative to the time averaged rate. 



SOME RESULTS ON TOTAL RATES 



For completeness and comparison we will briefly present our results on the total rates. Integrating 
the differential rates discussed in the previous section we obtain the total time averaged rate i?Oj 
the total modulated rate H and the relative modulation amplitude^, given by: 

R = Ro + Hcosa, R ^ Rq {1 + hcosa) (13) 

These are exhibited for zero threshold as functions of the WIMP mass in Figs \W\ and ?IU\ respectivelv. 
Some special results in the case of low WIMP mass are exhibited in Tables HlHIl In the case of non 
zero threshold one notices the strong dependence of the rime averaged rate on the WIMP mass. Also 
in this case the relative modulation h substantially increases, the difference between the maximum 
and the minimum can reach 20%. This however occurs at the expense of the number of counts, since 
both the time averaged and the time dependent part decrease, but the time averaged part decreases 
faster. So their ratio increases. This can be understood by noticing that the cancellation of the 
negative and positive parts in the differential modulated amplitide, see Fig. [1] becomes less effective 
in this case. 



(a) (b) 

a — !> 

Figure 10: The same as in Fig. IHlwith an energy transfer of 3 keVee (a) and 4 keVee (b). 




(a) (b) 

a — ^ 



Figure 11: The same as in Fig. |9]with an energy transfer of 5 keVee (a) and 6 keVee (b). 




Table I: Some the total event rates for some special WIMP masses and energy thresholds. The 
coherent nucleon cross section of cr„ = 10~^pb was employed. 



Eth 

(keVee) 


mwiMP 
GeV 


i?o(I) 

kg-y 


^(I) 
kg-y 


h(I) 


iio(Na) 
kg-y 


jy(Na) 
kg-y 


h(Na) 


i?o(NaI) 
kg-y 


i/(NaI) 
kg-y 


h(Nal) 





80 


16.3 


-0.311 


-0.019 


1.518 


0.028 


0.019 


14.0 


-0.259 


-0.018 





20 


25.8 


0.285 


0.011 


2.35 


0.050 


0.021 


22.2 


0.249 


0.019 





10 


18.4 


0.356 


0.019 


2.045 


0.046 


0.022 


15.9 


0.309 


0.019 


5 


80 


7.00 


-0.042 


-0.006 


1.133 


0.038 


0.034 


6.11 


-0.030 


-0.005 


5 


20 


2.72 


0.247 


0.091 


1.07 


0.065 


0.060 


2.47 


0.219 


0.089 


5 


10 


0.008 


0.001 


0.187 


0.303 


0.031 


0.103 


0.053 


0.006 


0.114 
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(a) (b) 

a — > 



Figure 13: The same as in Fig. [TUlwith a light target (Na or F). 




(a) (b) 
a — >• 

Figure 14: The same as in Fig. [TT]with a hght target (Na or F). 



DISCUSSION 

In the present paper we obtained resuhs on the differential event rates, both modulated and time 
averaged, focusing our attention on small energy transfers and relatively light WIMPS. We found 
that: 

• The relative modulation amplitude crucially depends on the WIMP mass. For small masses it 
exhibits normal behavior, but for large masses it changes sign (minimum in June). This effect 
is more pronounced in the case of heavy targets. 
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Q — ^keVee 

Figure 15: The same as in Fig. [7] for the target Nal. 
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(a) (b) 
a — >• 



Figure 17: The same as in Fig. [10] with a target of Nal. 
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(a) (b) 

a — ^ 



Figure 18: The same as in Fig. [TT]with a target of Nal. 



Table II: The same as in table U for cr„ = 2 x lO^'^pb relevant for the DAMA region. 





771 WIMP 


Ro{l) 


H(l) 


h(I) 


i?o(Na) 




h(Na) 


iio (Nal) 


//(Nal) 


h(Nal) 


(keVee) 


GeV 


kg-y 


kg-y 




kg-y 


kg-y 




kg-y 


kg-y 







80 


4.07 X 10* 


-776 


-0.019 


3.80 X 10^ 


70.2 


0.019 


3.50 X 10* 


-647 


-0.018 





20 


6.43 X lO'' 


712 


0.011 


5.87 X 10^ 


126 


0.021 


5.54 X 10* 


622 


0.011 





10 


4.61 X 10* 


891 


0.019 


5.11 X 10^ 


115 


0.022 


3.98 X 10* 


772 


0.019 


5 


80 


1.75 X 10* 


-105 


-0.006 


4.83 X 10^ 


95.0 


0.034 


1.53 X 10* 


-74.6 


-0.005 


5 


20 


6.80 X 10^ 


617 


0.091 


2.69 X 10^ 


162 


0.060 


6.17 X 10^ 


547 


0.089 


5 


10 


19.4 


3.62 


0.187 


757 


78.1 


0.103 


132 


15.0 


0.114 
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JTiwiMP — >GeV 

Figure 19: The total (time averaged) event rate in kg-y for a kg of target of ^^^I (a), of ^■^Na (b) 
and of Nal (c) assuming a coherent nucleon cross section cr,, = lO^^pb and a zero threshold energy. 



• The relative modulation amplitude depends somewhat on the energy transfer, especially at 
low transfers. 

• For WIMP masses less than 10 GeV, the difference between the maximum and the minimum 
could reach between 20% and 40% for a heavy target, but it is a bit less for a light target, 
depending on the enegy transfer. 

• The relative modulation amplitude for Nal is the weighted average of its two components, and 
in the low energy regime, between 1 and 6 keVee, it does not change much with the energy 
transfer. 

• Once it is established that one actually observes the modulation effect, the sign of the modu- 
lation may be exploited to infer the WIMP mass. 
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Figure 20: The total modulated event rate in kg-y for a kg of target of ^^^1 (a), of ^'^Na (b) and of 
Nal (c) assuming a coherent nucleon cross section = lO^^pb and a zero threshold energy. 



For low WIMP mass the total rates depend strongly on the threshold energy, especially for a heavy 
target. The relative modulation in the presence a threshold gets quite large {h « 0.2), but, unfortu- 
nately, this occurs at the expense of the number of counts. It is important to compare the relative 
total modulation in a least one light and one heavy target. For very low energy thresholds, if the 
signs are opposite, one may infer that the WIMP is heavy, jtiwimp > 100 GeV. 
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Figure 22: The total modulated event rate in kg-y for a kg of target of ^^^I (a), of ^■^Na (b) and of 
Nal (c) assuming a coherent nucleon cross section an — lO^^pb and a threshold energy of 5 keVee. 
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